Bimodality and hysteresis in systems driven by confined Levy flights 
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We demonstrate occurrence of bimodality and dynamical hysteresis in a system describing an 
overdamped quartic oscillator perturbed by additive white and asymmetric Levy noise. Investigated 
estimators of the stationary probability density profiles display not only a turnover from unimodal 
to bimodal character but also a change in a relative stability of stationary states that depends on 
the asymmetry parameter of the underlying noise term. When varying the asymmetry parameter 
cyclically, the system exhibits a hysteresis in the occupation of a chosen stationary state. 
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INTRODUCTION 

The Langevin description of an overdamped Brownian 
motion in a potential V{x) 



m = -v'{x) + at), 



(1) 



constitutes a basic paradigm to study the effects of fluctu- 
ations at the mesoscopic scale P, 0, Q . Here prime stands 
for differentiation over x and is a commonly assumed 
white Gaussian noise representing close-to-equilibrium 
fluctuations of the dynamic variable x. In contrast, in far- 
from-equilibrium situations, the Gaussianity of the noise 
term may be questionable, due to e.g. strong interaction 
with the surrounding "bath" 0, Q . A natural generaliza- 
tions to the Brownian motion are then 5-correlated Levy 
stable processes which can be interpreted as fluctuations 
resulting from strong collisions between the test parti- 
cle and the environment. In particular, the scale-free, 
self similar feature of Levy distributions gives rise to the 
occurrences of large increments of the position coordi- 
nates Ax during small time intervals causing a non-local 
character of the motion. Within the paper we address 
the problem of kinetics as described by Eq. ([1]) under 
the action of C(i) representing a stationary white Levy 
noise Q. Accordingly, the position of the Brownian 
particle subjected to additive white Levy noise is cal- 
culated by direct integration of E g. (fl ) with respect to 
the a-stable measure [1, 0, H, H, flQl ll| La,i3{s), i.e., 
x(t + At) = x(t) -V {x(t))At + {At)'^/°'C, where C is dis- 
tributed according to the a-stable Levy type distribution 
ia,/3(C) M — 0) whose representation [y, H, [l3| is given 
by the characteristic function defined in the Fourier 

space 0(fc) = J^(La,;3(C;o-,M)) = '^Ce'''^ia,/3(C; o", 
for a 1 



(j){k) = exp -cr"|fc|" - i^sgn(fc) 

and for a = 1 
4>{k) — exp 



tan J + i^k 



-a\k\ l + i(3-sgn{k)ln\k\ 

TT 



(2) 
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The stability index a, determining tails of the probability 
density function (PDF) takes values a S (0, 2], the skew- 
ness of the distribution is modeled by the asymmetry pa- 
rameter /3 £ [—1, 1]. Indices a and /3 classify the type of 
stable distributions up to translations and dilations 6, 8]. 
Two other parameters of scaling a G (0, oo) and location 
/i G (— cx), cx)) can vary, although replacing C ~ M ^^'^ '^C 
in the original coordinates, shifts the origin and rescales 
the abscissa without altering function La, /){(,) ■ For sim- 
plicity, we will restrict here to strictly stable Levy noises 
with ^ = [i, [i]. Generally, for /3 = = PDFs are 
symmetric while for (3 — ±\ and a G (0, 1) they are to- 
tally skewed, i.e., C, is always positive or negative only, 
depending on the sign of asymmetry parameter (3 (cf. 
Fig. [2|). Asymptotically, for C ^ oo with a < 2, stable 
PDFs behave as p{C,) oc causing divergence of 

moments {Q') = oo for u > a. The asymmetry is re- 
flected in a biased distribution lim^^oo p(|zj>'i;) ~ ^~2^- 
Equivalent to the stochastic ordinary differential 
Eg . IT]) is a fractional Fokker-Planck equation 



15l . Il6l (FFPE) for the probability distribution func- 



tion 
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'|fc|"+i/3CT"/fc|fc|' 



tan ■ 



-A)"/2p(x,t) 
TTa d 



i-A) 



(q-1)/2 



p{x, t) (4) 



with p{x, t) representing the probability density func- 
tions for finding a particle at time t in the vicinity of 



X, (j){k,t) = {e'Kp[ikx{t)]) = ^exp ik C(s)(is J stand- 
ing for the characteristic function of the stable pro- 
cess and —A"f{x) denoting fractional Laplacian [Tsl. [l9l| 

-(A)"/2/(a;) = (\k\°'f{k)^, with a = 2 correspond- 
ing to the standard Brownian diffusion case. The addi- 
tion of the potential force —V'{x) to Eq. ^ adds the 
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classical drift term \V' {x)p{x,t)] to Eq. dH). In the 
approach presented herein, instead of solving Eq. ([4]) , in- 
formation on the system is drawn from the statistics of 
numerically @,[1,[20] generated trajectories satisfying the 
generalized Langevin equation ([T]). At a single trajectory 
level sampled from the stochastic dynamic study of the 
problem, the initial condition for Eq. ([1]) has been set 

to x(0) = W[— 1,1], i.e., initial position of the particle 
is drawn from the uniform distribution over the interval 
[—1,1]. For simplicity. E g. (fll) has been studied in di- 
mensionless variables 2l|, |22| with additionally setting 
a — 1. The time independent potential V{x) is assumed 
to be of the form V{x) = a;'*/4 which guarantees the con- 
finement of the trajectory x{t) within the potential well 
[21, 23, 24] leading to a finite variance of the stationary 
PDF. For the general a-stable driving and the quartic 
potential the stationary PDF fulfills 



d^P{k) 



= sgnfclfcl"-ip(fc) 



z/3tan^lfc|"-ip(fc) (5) 



where P{k) stands for a Fourier transform of the station- 
ary PDF, P{x) = \\Ta.t^oo P{x,t). Analytical solutions 
of Eq. © can be readily obtained for a Gaussian case 
V4) and for a Cauchy additive 



a = 2: P{x) oc exp( 
noise a = 1: P{x) ■ 



1/[^(1 - a;^ + x^)] 12l|, lai, [241. 
They display a perfect agreement with profiles of PDFs 
obtained by numerical simulations of E g. (f Tl) performed 
with the Janicki-Weron algorithm [^, [^, [SsTlia] . 



RESULTS 

Numerical results were constructed for a time step of 
integration At = 10"'^, simulation length T^ax — 10 with 
an overall statistics of = 5 x 10'' realizations. To check 
whether results are influenced by the length of simula- 
tions, results for various T^ax (Tmax = 10, Tmax = 15) 
were compared showing consistency of estimated PDFs 
for both values. Further details on values of parame- 
ters are included in the text underlying the figures. Nu- 
merical examination of Eq. ([1]) allows for construction of 
PDF estimators for the whole range of parameters a and 
/?. Likewise, by direct integration of Eq. ^ it is also 
possible to investigate time- dependent PDFs [l^] and 
noise-induced bimodality of the probability distribution 
23,|2i|,[2l|23. 

The change of a (for /3 / 0) from values greater than 1 
to values smaller than 1 results in change of location of a 
modal value of stable densities, see left vs. right panel of 
Fig. [2l The shift of modal values is reflected in properties 
of stationary states. For a > 1 with /3 < the modal 
value is located for a; > (right panel of Fig. [3|) while 
for a < 1 with < it shifts to the negative interval 
X < (left panel of Fig. [31). For /3 < modal values are 
located on the opposite side of the origin than for /3 > 0. 
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FIG. 1: Stationary probability density functions (PDFs) for 
a = 2 (Gaussian case, unimodal distribution) and a = 1.0 
(Cauchy case, bimodal distribution) with corresponding an- 
alytical solutions. Numerical results were constructed for 
At = 10"^, Tmax = 10 and averaged over A = 5 x 10* re- 
alizations. 
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FIG. 2: Sample a-stable PDFs with a = 0.9 (left panel) and 
Of = 1.1 (right panel). For P = Q distributions are symmetric 
and become asymmetric for /3 = ±1. The support of the 
densities for the fully asymmetric cases with /3 = ±1 and 
a < 1 (left panel) assumes only negative values for (3 = —1 
and only positive values for 13=1. Note the differences in the 
positions of the maxima for a < 1 and q > 1. 



Therefore, the position of the modal value can be moved 
from the positive to negative real lines by change of a 
while (3 is kept constant (left vs. right panel of Fig. [3]) 
or by change of (3 to —(3 with a preset value of a (cf. 
different rows in Fig. [3]). 

The fact that changes in (3 can change the location of 
the modal value suggests that periodic changes in this pa- 
rameter can lead to a phenomenon resembling dynamical 
hysteresis 27, 28, [2^. In order to register a hysteretic 
behavior of the system, we have performed an analysis of 
trajectories of Eq. ([T]) based on a two-state approxima- 
tion. For that purpose, we have defined an occupation 



3 



1 r- 
3 - 

2 - 
1 - 

L 
-3 



^ 3 
& 2 
°- 1 






0.6 
0.4 
0.2 


0.6 
0.4 
0.2 




„ 0.6 r— 
& 0.4 - 
°- 0.2 - 

L 

3 -3 

„ 0.6 ^ 
& 0.4 - 
°- 0.2 - 


3 -3 



FIG. 3: Stationary probability distributions p{x) for a = 0.9 
(left column) and q = 1.1 (right column). Various rows corre- 
spond to the various values of starting from the top panel: 
f3 ^ I (top row), /3 = 0.5, f3 = 0.01, (3 = (bottom row). 
Due to symmetry, results for negative (3 can be obtained by 
the reflection of results for /3 ^ along a:: = line. Results 
were constructed for At = 10~*, Tmax = 10 and averaged over 
N = 5 X 10** realizations. 



probability of being in left /right state according to 
pt(left) = prob{x(i) < 0} = / p{x,t)dx 

J —QO 



= 1 -_Pi (right). 



(6) 



Transition between the states is induced by time depen- 
dent asymmetry parameter (3 which is periodically mod- 



ulated over time, i.e., f3 = /Jmaxcosrit = /3, 



max cos rp^_^ 



In 



the Fig. m values of p(left) for various To (with Pniax = 1) 
are presented. Stable random variables for a < 1 with 
|/3| = 1 take only positive/negative values, depending on 
the sign of f3. Therefore higher level of saturation is ob- 
served for a = 0.9 (left panel) than for a = 1.1 (right 
panel), i.e., when = 1 probability of being in the 
left/right state for a = 0.9 is higher than for a = 1.1 
(right panel). The direction of the hysteresis loop is a 
direct consequence of the fact that changes in /3 move 
modal values from positive/negative real line to the neg- 
ative/positive real line. Due to the initial condition im- 
posed on x(0), the starting point for each hysteresis loop 
is (0,0.5) and the first part of the curve describes ap- 
proaching to the proper hysteresis loop. With decreas- 
ing O (increasing Tq) area of hysteresis loop decreases 
because the system has more time for relaxation and re- 
sponse to the changes of /3. For large Q (small Tq) the 
response of the system is more delayed and consequently 
area of the hysteresis loop is larger. 

Finally the influence of /9max on the shape of hystere- 
sis loop has been examined. In Fig. [5] hysteresis loops 
for various /^max (with Tji = 9) arc presented. Smaller 
/3max makes stable distribution less skewed and as a con- 
sequence, less probability mass becomes located in the 
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FIG. 4: (Color online) Hysteresis loops for a = 0.9 (left panel) 
and Q = 1.1 (right panel). Different lines correspond to the 
different values of Tn, amplitude of /3 /3max = 1.0. Results 
were constructed for At = 10~^ and averaged over 5 x 10* 
realizations. Various loops corresponds to different values of 
the driving period Tn- Tn = 35, Tn = <d and Tn = 1.5 (from 
the inside to outside). 



left state and the effect of saturation is less visible. Fur- 
thermore, with decreasing /^maxj loops become more oval 
and finally for /Jmax — the hysteresis phenomenon dis- 
appears. 
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FIG. 5: (Color online) Hysteresis loops for a = 0.9 (left panel) 
and Q = 1.1 (right panel). Different lines correspond to the 
different values of /3max, driving period Tn = 9. Results were 
constructed for At = 10~^ and averaged over 5 x 10* realiza- 
tions. Various loops corresponds to different values of /Jmax: 
/3max = 1.0, /3max = 0.5, /3max = 0.3 and /Jmax = 0.1 (from the 



top to bottom). 



SUMMARY AND CONCLUSIONS 

The modulation of stable noise parameters modify 
shape of stationary densities corresponding to Eq. ([T]) . In 
particular skewed noise characterized by nonzero asym- 
metry parameter (3 can induce asymmetry of station- 
ary states in symmetric potentials. Furthermore, totally 
skewed stable noises (a < 1 with = 1) could move 
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the probability mass to one side of the x-axis making 
stationary states totaUy skewed. In such situations, the 
whole probability mass can be located on the left hand 
side or on the right hand side of the origin x = de- 
pending on the sign of the asymmetry parameter /3. For 
a > I with a nonzero /3, stable noises still are asymmet- 
ric. Nevertheless, even |/3| = 1 is not sufficient to induce 
totally skewed stationary states. Consequently, dynam- 
ical hysteresis loops induced by cyclic variation of f3 for 
a < 1 and a > 1 are characterized by various level of 
saturation, see left vs. right panel of Fig. 01 

The stationary densities depicted in Fig. [3] are recorded 
in the system described by Eq. ^ which is subjected 
to the action of stable noises. The very same sta- 
tionary densities can be also observed in the equilib- 
rium system perturbed by the white Gaussian noise, i.e., 
= -Kff(3^'0 where Vesix,t) = -log[p(x,i)] 

and {£.{t)£,{t')) = S{t — t'). The dynamical hysteresis de- 
tected in the model described by Eq. ^ emerges as a 
consequence of periodical modulation of the asymmetry 
parameter (3. In the effective potential model, periodical 
modulation of the asymmetry parameter corresponds to 
the periodical modulation [30] of the effective potential 
Ves{x). Nevertheless both models are not fully equiv- 
alent. They have the same one-point densities p{x,t), 
while other characteristic of the process {x{t)} are dis- 
tinct [22]. Therefore, examination of stationary densities 
itself is not a fully conclusive method for discrimination 
of underlying types of noises. 

The dynamical hysteresis can be also observed in 
the generic double well potential model subject to the 
joint action of the deterministic periodic modulation and 
stochastic a-stable fluctuations 20]. In such a case, how- 
ever, the shape of the dynamical hysteresis loop is af- 
fected both by the character of the noise and by the type 
of periodic modulation 20] . 
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